Equation of state of a weakly interacting two-dimensional Bose gas studied at zero 
temperature by means of quantum Monte Carlo methods 
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The equation of state of a weakly interacting two-dimensional Bose gas is studied at zero temper- 
ature by means of quantum Monte Carlo methods. Going down to as low densities as no? oc 10~ 100 
permits us for the first time to obtain agreement on beyond mean-field level between predictions of 
perturbative methods and direct many-body numerical simulation, thus providing an answer to the 
fundamental question of the equation of state of a two-dimensional dilute Bose gas in the universal 
regime (i.e. entirely described by the gas parameter na 2 ). We also show that the measure of the 
frequency of a breathing collective oscillation in a trap at very low densities can be used to test the 
universal equation of state of a two-dimensional Bose gas. 



The calculation of properties of weakly-interacting 
quantum gases was historically a very important and 
productive task as it led to the development of pertur- 
bative methods, such as the Bogoliubov diagonalization 
approach, Fcynman diagrams, etc. The equation of state 
of a three-dimensional Bose gas dates back to 1957 when 
the beyond mean-field correction was obtained by Lee, 
Huang, and Yang [lj. In the beginning of 1960s the 
equation of state of a dilute one-dimensional Bose gas 
was obtained by Girardeau [2|. From this perspective it 
is rather surprising that even until now the equation of 
state of a dilute two-dimensional (2D) gas remains an 
open question. In the last thirty years, theproblem has 
been addressed in different studies 0, 0, H, 0, 0, @, @, 
11 , 12l 13L 14, li| but the results obtained are often in- 
complete and even contradictory, which indicates that 
the two-dimensional problem is extremely complicated. 
Indeed, a peculiarity of two-dimensional systems is that 
the coupling constant depends on the density, contrarily 
to three- and one- dimensional cases, where it is entirely 
defined by the s-wave scattering length a and the particle 
mass to. This feature makes the analytical derivation of 
the equation of state much more involved. 

When perturbative theories fail one can resort to ah 
initio numerical simulations. The correctness of pertur- 
bative equations of g round state has been checked in 
three dimensions [161 ] and in one dimension (see, e.g. 
Ref. [nl) using diffusion Monte Carlo (DMC) methods. 
Analogous attempts ID, lj| to find numerical agreement 
with beyond mean-field terms in 2D systems have not 
succeeded as in this case the expansion parameter has 
logarithmic dependence on the gas parameter na 2 . In 
particular, densities as low as na 2 oc 10 -5 [20j, 10 -7 [HI 
and 10 -6 [n| have been reached. Our estimation below 
shows that for confident testing of the beyond mean-field 
(BMF) terms the density should be smaller than 10 -69 , 
which obviously is an extremely challenging task. In this 
work we do an effort towards calculations at considerably 
smaller densities in order to test different expressions for 



beyond mean-field terms present in the literature. 

We use the diffusion Monte Carlo (DMC) method 
to address the problem. DMC solves stochastically 
the Schrodinger equation providing exact results for the 
ground state energy within controllable statistical errors. 
The construction of the guiding wave function is as in 
Ref. 2l| • In Table U we report the energy per particle 
in the thermodynamic limit as a function of the density. 
The scattering at low densities is universal: it is indepen- 
dent of a particular choice of interaction potential, and 
is described by a single parameter, the s-wave scattering 
length. Although the calculations are done for the dipo- 
lar interaction there is no difference between results for 
different interaction potentials for densities nr 2 <C 10~ . 
For convenience we keep the name of "dipoles" in order 
to distinguish between different series of calculations. 

The leading contribution in the 2D equation of state 
is described by mean-field Gross-Pitaevskii (GP) theory 
in which all particles are considered to be in the con- 
densate. The chemical potential is proportional to the 
density /1 = g2r>n and the coefficient of proportionality 
is the coupling constant. Unlike three- and one- dimen- 
sional systems, where the coupling constant is indepen- 
dent of the density, in a two-dimensional system 52D itself 
depends on density making the construction of a precise 
theory very complicated. The leading contribution to the 
purely 2D coupling constant was first obtained by Schick 
|3| in 1971, who made use of the Beliaev method 22]. 
Contrary to three- and one- dimensional cases, the de- 
pendence of the coupling constant on the s-wave scatter- 
ing length a is very weak. Indeed, a enters under the 
logarithm, 
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Recently, it has been shown that a rigorous derivation of 
a two-dimensional Gross-Pitaevskii theory leads to the 
same result fl2|. 

Before discussing beyond mean-field terms it is nec- 
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Table I: DMC energy per particle for dipolar interaction as 
a function of the density nr 2 . obtained by extrapolating the 
results for N = 100 — 800 particles to the thermodynamic 
limit. Conversion to units of s-wave scattering length a can 
be easily done according to the relation a — 3.17222ro. The 
energy is given in units of h 2 jrar\. Statistical errors on the 
last digit are shown in parentheses. 



essary to note that the usual mean-field relation be- 
tween energy and chemical potential is accurate only in 
the leading term, while in three- and one- dimensional 
systems the MF relation is exactly linear. Indeed, the 
chemical potential is proportional to the coupling con- 
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which is found by integrating the chemical potential, is 
E MF /N = 47rr(0,2|lnna 2 |)?i 2 /(mna 4 ) where F(a,x) is 
the incomplete gamma function. The energy expansion 
in the dilute regime na 2 — ► can be obtained from the 
large argument expansion of r(0,x) as follows: 
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We see that there are differences between Eq. ([2]) and the 
usual MF expression E MF /N — gn/2. The differences 
contain logarithmically small terms that exceed the ac- 
curacy of the mean-field theory. However, in the study 
of beyond MF effects terms of this order are significant. 

The dominant beyond mean-field terms were obtained 
by PopovQ in 1972 (see also his book [23[). He obtained 
a recursive expression relating chemical potential \x and 



density n. At zero temperature his expression reduces to 

(3) 
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where Eq is of the order of ti 2 /mrg, r$ being the range of 
the interaction potential. We write the last relation in- 
troducing an unknown coefficient of proportionality C\: 
c = Cih 2 /ma 2 . By solving Eq. iteratively one ob- 
tains the following expression for the chemical potential 
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Mathematically, the leading beyond MF term in ijl and 
E/N is a double logarithm ln|lnna 2 |, and if the den- 
sity is low enough inclusion of this term should be suffi- 
cient to reproduce the energy correctly. From a practical 
point of view this term works at such extremely low den- 
sities; that inclusion of the double logarithm term alone 
might even lead to larger deviations in the energy. In- 
deed, the second subleading term is some constant of the 
order of ln47r. To estimate a typical density of applica- 
bility of a double-logarithm term we ask that the term 
ln | ln na 2 1 be m times larger than typical constant term 
ln47r. This leads to an estimation of a characteristic den- 
sity na 2 w e~( 47r ) . Asking for a factor m — 2 (i.e. 50% 
accuracy) we get an extremely low density na 2 « 10 -69 , 
asking for m = 3 (i.e. 33% accuracy) the density drops to 
a quite unrealistic number na 2 ~ 10~ 862 . We note that 
the universal regime where the energy depends only on 
the gas parameter na 2 extends to much larger densities, 
such as na 2 « 10~ 6 . This situation should be contrasted 
with three- or one-dimensional cases, where the univer- 
sal regime is perfectly described by the first beyond MF 
correction. We note that problems due to a logarithmic 
correction arise also in four dimensions 24]. Another 



important peculiarity of two dimensional world is that 
second and third beyond MF corrections have different 
signs. Furthermore, it follows that such corrections can 
compensate each other at some relatively "large" density. 
In fact, this effect has been observed in a full many-body 
calculation (l9j at density na 2 w lx 10~ 6 . Notwithstand- 
ing the energy in this point is reproduced correctly by 
the mean-field theory, other properties (for example, the 
condensate fraction) are not described precisely. 

The complexity of the problem explains the large num- 
ber of works dedicated to the analytical study of the 
equation of state of dilute 2D Bose gas. A summary 
of these works is presented in Table [U where we report 
the expression for the energy per particle. Almost all be- 
yond MF theories agree on the presence of the ln | ln na 2 1 
term but differ in the constant, or explicitly discuss that 
the constant term can not be obtained within the used 
approach. Previous diffusion Monte Carlo calculations 
la . Il9j ] have shown the existence of a universal equation 
of state, but even at the smallest considered densities 
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Table II: Literature overview: equation of state of a dilute 
two-dimensional Bose gas. Fifth column, terms of the low- 
density expansion appearing as straight lines in Fig. [T] in the 
limit of large In | In na 2 \ . 
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Figure 1: (Color online) Non-universal beyond MF corrections 
2Tvh 2 n/[mE/N}- \ In na 2 1 in the energy per particle as a func- 
tion of the double logarithm of the gas parameter In | In na 2 \ . 
Symbols, DMC results: up triangles, hard disks[l8j], down tri- 
angles dipoles. Lines, different equations of state, see TablelHl 
Black dashed line, Eq. JB). 



none of the equations of state given in Table |TT] was able 
to reproduce the energy correctly. 

In the present DMC calculation at ultra-low densities 
we reach the regime where the analytical expansions from 
Table [U are applicable and therefore it is possible to test 
them. According to the majority of theories the complete 
expression for the energy can be generalized as E/N — 
2ith 2 n/[m{\ lnna 2 + In | lnna 2 — const + ...)]. Then, the 
beyond MF terms should be tested by inverting the en- 
ergy and subtracting the MF term, according to the com- 
bination 2Trh 2 n/[mE /N] — | lnna 2 | = In | lnna 2 | — const. 
Figure Q] shows DMC results of this combination as a 
function of the double logarithm of the gas parameter 
In | In na 2 1 for hard disk and dipolar interaction poten- 
tials. The analytical results from Table HT1 are drawn with 
lines. The majority of approaches predict the first be- 
yond mean-field term as In | In na 2 1 that would correspond 
to a linear behavior with unitary slope in Fig.[TJ The only 
two exceptions are the theory of Hines et al. [6| and a fit 
to DMC results of hard disks [18| . The region of applica- 
bility of the latter fit is limited to 1.5 < In | lnna 2 | < 2.8. 
Our results for dipolar interaction potential are compat- 
ible with the linear dependence and unitary slope (for 
the largest x = In | In na 2 | > 4 the numerical simulations 
are very difficult and the error AE in the estimation of 
the energy is exponentially amplified in the quantity of 
interest as (nAE)/E 2 w e x {AE)/E). The third beyond 
mean-field correction corresponds to a constant shift in 
Fig. [TJ Here there are the major discrepancies between 
theories. The reason for that is the complexity of per- 
turbation theory in two-dimensions. One has to evaluate 
correctly all the terms contributing to a given level of 
accuracy. Just to name some sources for this constant, 
already from the recursive expression ^ one sees that 
the self-consistent evaluation of /x contributes to its value. 



Also the effect of the quantum fluctuations (depletion of 
the condensate density), the difference between cut-off 
length, the range of the potential and the s-wave scatter- 
ing length contribute on this level. The Euler's gamma 
constant 7 appearing in some of the theoretical results 
is related to the dependence in the ladder approximation 
of the chemical potential on the scattering amplitude [22|. 
which in turn contains 7 in its short-range expansion^]. 

Our DMC results are in agreement within error bars 
with the constant shift — ln7r — 1/2 — 2j obtained in 
Refs. [Tl|, El, EH for ln|lnna 2 | > 3. Cherny and 



Shanenko do expansion in terms of a dimcnsionless in- 
medium scattering amplitude u obeying the equation 
l/u + In it = — ln(7rna 2 ) — 27. The dimensionless en- 
ergy e = Em/ (2irh 2 nN) is then expanded as a se- 
ries e = u + u 2 /2 — u 3 + ... We test the accuracy of 
this description and find that the agreement with the 
first two terms of this expansion is notably good. The 
deviations enter on the level of the third beyond MF 
term. We apply a fitting procedure to find the coeffi- 
cient in front of u 3 with the % 2 criterion; the result is 
mE/(2Trh 2 Nn) = u + u 2 /2 - 2.0(l)u 3 . This fit describes 
correctly the energy in the universal regime (where re- 
sults for dipoles and the hard disks coincide) and at den- 
sities In I In na 2 \ > 2 for hard disks. We note that the 
hard disk interaction potential is completely described 
by only one parameter and its equation of state is ex- 
pected to be the most universal. 

In terms of the gas parameter our best perturbative 
description for the chemical potential is given by 
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Figure 2: (Color online) Main plot: energy in units of 
MF energy E™ F /N = 2tt h 2 n/\m In na 2 \ as a function of 
the gas parameter, inset : square of the lowest breath- 
ing mode frequency, Q 2 as a function of the LDA param- 
eter N 1 ^ 2 r 2 ) /a 2 lo . Thin dashed line (inset), mean-field lin- 
ear equation of state ^ ¥ oc n, thin solid line (inset and 
main plot): mean-field expression of Schick for the energy 
E^ F , thick solid line (main plot): mean- field expression 
Ef F /N = 47rr(0, 2| In na 2 \)h 2 / (rnna 4 ) dash line with one dot: 
first beyond MF correction, dash line with two dots, second 
BMF correction, short-dashed line: three BFM terms, circles: 
DMC results for dipoles, squares: DMC results for hard disks. 



where Cf = — In ir — 27 — 1 = —3.30... and Cf = — hi7r — 
27 + 2.0(1) = —0.3(1). The corresponding expression for 
the "universal" energy per particle is 

E 2-Kh 2 n/m 

N ~ M 21 1 1 II 21 1 nE , ln|lnna 2 |+C 2 E . * ' 

1 v I In na 1 \ + In | In na 1 \ + Cf H | ln wa2 — *- + . . . 

where Cf = Cf + 1/2 and Cf = Cf + 1/4 (in the mean- 
field case, described by Eqs. ([TJ) and relation between 
Cf and Cf is similar, while relation between Cf and 
Cf is different due to double-logarithm term present in 
Eq. ©). 

We illustrate the convergence of the series in Fig. [2] 
where we plot the quotient between the energy and mean- 
field prediction of Schick[3]. One sees that contrary to 
three- and one- dimensional cases, the mean-field pre- 
diction is not precise even at extremely dilute values of 
the gas parameter [l9jj. Moreover, differences between the 
expression of Schick and the integrated Schick's chemical 
potential matter at the considered densities. Inclusion 
of the leading beyond mean-field correction (i.e. double- 
logarithm term ln | lnna 2 |) does not improve the descrip- 
tion for densities na 2 > 10~ 100 . This term is negative 
and it leads to an underestimated energy. The second 
BMF correction (Cf and Cf terms in Eqs. P5]l ) has a 
different sign and effectively cancels the first BMF term 
at the density na 2 ~ 10 -6 . Inclusion of this term permits 



to recover the DMC results up to na 2 « 10~ 30 . Inclusion 
of the third BMF term, Eq. ([6]), extends the region of ap- 
plicability of the series expansion up to relatively "high" 
densities na 2 « 10~ 3 . 

The correctness of the series expansion of the 2D equa- 
tion of state can be verified in high precision experiments 
measuring the frequency Q of the lowest breathing mode. 
We note that the presence of a tight transverse harmonic 
confinement modifies the coupling constant reducing the 
effectively s-wave scattering length a in (TT]) by an ex- 
ponentially small factor ~ exp{— aho/a3D}, where 
ciho is oscillator length and 030 is the s-wave scatter- 
ing length [26]. This makes ultradilute 2D densities be 
experimentally accessible exploiting confinement induced 
resonance. We use local density approximation (LDA) to 
describe a trapped gas and show our predictions for £1 in 
the inset of Fig. [2j As one can see, the effect is large 
and therefore the observation of f2 provides a sensitive 
experimental tool for testing the equation of state. 

To conclude, we study the zero temperature equation 
of state of weakly interacting 2D Bose gas down to ultra 
dilute densities na 2 oc 10 -100 and confront the obtained 
energy with a large number of different analytical ex- 
pressions for the energy. All theories agree at mean-field 
level and the majority agrees on the first beyond MF 
term, which, still does not provide an accurate quanti- 
tative description even at such low densities. A good 
agreement between perturbative series and results of the 
present ab initio calculations is found for the first time. 
The third beyond MF term is fitted and the correspond- 
ing equation of state is explicitly written in terms of the 
gas parameter. Finally, we suggest future experimental 
work to measure the frequency of a breathing collective 
oscillation in a trap at very low densities since our re- 
sults prove that this frequency is a very sensitive tool for 
verifying the universal 2D equation of state. 
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